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^ ■ 

^ . Abstract 

We consider Wick's Theorem for finite temperature and finite volume systems. Work- 
1^ ■ ing at an operator level with a path ordered approach, we show that contrary to claims 
in the literature, expectation values of normal ordered products can be chosen to be zero 
^ , and that results obtained are independent of volume. Thus the path integral and opera- 
te ' tor approaches to finite temperature and finite volume quantum field theories are indeed 
seen to be identical. The conditions under which normal ordered products have simple 
^ symmetry properties are also considered. 

o ■ 

^ ■ 1 Introduction 

Qh! Early developments of thermal field theory were based on a path ordered operator ap- 
^ I proach |]T| , and Wick's theorem played a crucial role in the development of its perturbation 
theory [ffl, 0]. More recently, other approaches to the subject have been developed - in 



X 



particular, path ordered methods have been extended, often by using a path integral ap- 
proach 1^, ^, ^, , and Thermo Field Dynamics (TFD) in which a doubhng of the states 
of the system is introduced 0, H, ||, |3, ||, 3 . It is important to understand to what degree 
these various approaches are equivalent [^, 

However, a closer inspection of the original literature and subsequent standard texts 
turns out to be rather confusing regarding the issue of Wick's theorem and the expectation 
value of normal ordered products. In particular there are claims that there are finite 
volume corrections to the standard thermal perturbation theory. If this was true then 
in some important cases, such as in the study of the quark-gluon plasma formed in 
relativistic heavy ion collision experiments, there would be important new corrections 
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to some existing calculations. In this paper we concentrate on path-ordered operator 
approaches and we clarify the confusion found in the literature. These results have a 
direct translation into a path integral language which we consider in the penultimate 
section. In terms of TFD the principles are similar but the implementation is rather 
different. However, the case of TFD is considerably easier because of its similarity with 
zero temperature field theory and the literature reflects this. We will consider TFD in 
our conclusions. 

Central to our discussion is Wick's theorem which we refer to as an operator iden- 
tity rather than a relation between expectation values. It relates time ordered products 
of fields to normal ordered products and contractions. At zero temperature, the use 
of Wick's theorem revolves around the fact that normal ordered products are defined 
such that annihilation operators are always placed to the right of creation operators. 
This ensures that vacuum expectation values of normal ordered products, A^[. . .], vanish: 
(0|iV[aa'^] |0) = 0. However with the ensemble averages used in thermal field theory, prod- 
ucts of fields normal ordered in the same way are no longer zero: 

Since in calculations one considers Green functions, expectation values of time-ordered 
functions, it is useful but not crucial to show that thermal expectation values of normal 
ordered products are zero: all that is actually required is a relation between Green func- 
tions of different orders. It is this 'expectation value version' of Wick's theorem which is 
proved by Thouless [|l^ for real relativistic scalar fields. Using rather different techniques, 
Gaudin proved the same 'expectation value version' in a rather more transparent way 
using the cyclicity of the thermal trace. This is the approach also used in P, [12], |l3l; it 
avoids the need to discuss normal ordering at all. Volume does not enter into the argu- 
ment either. In more modern texts 0, ^ [T^ which rely on path integral methods, the 
expectation value version of Wick's theorem is inherent in the notation - it is essentially 
encoded as Schwinger-Dyson equations for the free field case. 

We now return to the issue of normal ordering which was not properly dealt with in 
the original paper of Matsubara. The area has been studied in subsequent literature, 
and seems to be a source of confusion. By using a definition of normal ordering different 
from that of standard field theory, Thouless [|I^ proves that thermal expectation values 
of normal ordered products of non-relativistic fermions can be set to be zero. Fetter and 



Walecka |12| use Gaudin's argument and claim that "Unfortunately, . . . the ensemble 



average of normal ordered products is only zero at zero temperature" contrary to |jTO|. 
Gross et al. |]T3| use the same approach. 

The situation is further confused by various authors' considerations of finite volume 
corrections. Matsubara suggests that the expectation values of normal ordered terms 
may be zero only in the infinite volume limit. This argument is given in greater detail 
in Abrikosov et al. and Lifschitz and Pitaevskii fl^. If this was true, there could be 
important corrections to standard perturbation theory, coming from the usually neglected 
normal ordered terms, when working with finite volumes systems. As the quark-gluon 
plasma formed in relativistic heavy ion collisions fills a comparatively small volume, such 
finite volume corrections might be of great importance in practical calculations. Yet where 
are such volume terms hiding in the derivations using the properties of the thermal trace 



11, 12, 13 1 and in path integral derivations p, a, H? 



The finite volume argument is worth repeating here for reference later on. For sim- 
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plicity, consider a non-interacting equilibrium system of relativistic neutral bosons at 
temperature T lying in a volume V. Let iplx) be the hermitian field operators in the 
Heisenberg picture. We wish to calculate Green functions, and as a specific example 
consider the four-point thermal Wightman function 

>Vl,2,3,4 =< ^iXi)^{x2)ijiXs)^ix4) > • (1.1) 

The angular brackets denote a thermal expectation value (which is defined in the next 
section), and Xi are different space-time points. In momentum space, Wi, 2,3,4 contains 
terms of the form ^ 

772 imim ^ 4i«L«k3ak4 > exp(3^), (^-2) 

ki k2 ks k4 

where and cik are creation and annihilation operators, the kj's are discrete momenta, 
and 3^ is a sum of scalar products of 4 momenta and space-time coordinates. In the 
literature it is claimed that there are three separate contributions to ( |1.2| ) [IB], . The 



first 2 contributions arise when the momenta are equal in pairs, that is when ki = ks, 
k2 = k4 and ki = k4, k2 = ks. For example, in the former case ( |1.2| ) becomes 

^ < 4i42«kiak2 > exp(3^). (1.3) 

^ ki k2 

The final contribution arises when all the momenta are equal: 
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< 44«kak > exp(3^). (1.4) 



According to some standard texts (for example 0, |T7|); "expression ( [1.3|) has a special 
feature which distinguishes it from ( p..4|) , i.e. the number of factors in 1/V is the same 



as the number of summations, whereas in ( |1.4| ) this number of factors is larger". Thus 



( |1.4D has "an extra factor of 1/V and will vanish as ^ ^ 00". These texts conclude that, 
once the infinite volume limit has been taken, "we can actually average over pairs of 
operators" . The implication is that only in this limit does the four-point Green function 
(or in general an ra-point one) split into a sum of products of two-point Green functions. 

In this paper we prove that the above argument is incorrect and that any n-point 
Green function can be decomposed exactly into two-point Green functions, whatever V. 
We also prove explicitly that thermal expectation values of normal ordered products of 
any type of field may always set to zero, contrary to claims elsewhere [|12|. The symmetry 
properties of normal ordered products are considered. 

The paper is organised as follows. In section 2 we set up our notation and definitions. 
This is essential as we will be discussing different types of normal ordering and we wish 
to track volume factors very carefully. In section 3 we split fields into two parts at an 
operator level and define normal ordering and contractions appropriately. We choose a 
split which guarantees that the thermal expectation value of normal ordered products 
of any two field operators is always zero and hence that the normal ordered products 
have simple symmetry properties. We then prove in section 4 that this is sufficient to 
ensure that the thermal expectation value of any four-point normal ordered product also 
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vanishes, independent of volume. With Wick's theorem, this then shows directly that the 
four-point Green function decomposes exactly into two-point Green functions. In section 
5 we generalise this result from four-point to ra-point Green functions by using more 
powerful methods. One is the usual path integral approach to thermal field theory. The 
other is based on the cyclic property of the trace which defines the thermal expectation 



value |TT|, |12|. We use these, with Wick's theorem, to show that the thermal expectation 
value of all n-point normal ordered products vanish. We also note the symmetry of the 
n-point normal ordered products. 



2 Definitions and Notation 



For simplicity we consider only free field theories and work in the Heisenberg picture. 
Our results are completely equivalent to working with an interacting theory in the in- 
teraction picture. There the operators evolve according to the free Hamiltonian and the 
interactions are encoded in the S matrix. When expanded, the S matrix leads one to 
consider time-ordered products which are of the same form as those considered here. At 
finite temperature, one has to be rather more careful about the use of an 5* matrix, but 
one can still follow a similar algorithm so that our results can be applied to practical 
calculations at finite temperature. 

Field operators, tp{x), are given in terms of particle and anti-particle creation and 
annihilation operators, aj^, dk and foj^, 6^ respectively; 



OkC 



-ik 



" + bl 



Ak.x 



and ujI_ 



+ for relativistic field theories. 



(2.1) 



(2.2) 



For 



The metric is diag(-|-l, —1, —1, —1) 
non- relativistic ones tUk = k^/2m, and &k = 0. A real relativistic scalar theory has 
Ok = For fermionic theories, ay_ and 6k implicitly carry spin indices, and there is an 
implicit sum over spins as well as over momenta. Annihilation and creation operators 
must also be multiplied by the correct solutions to the Dirac equation if the theory is 
relativistic but this does not effect our results. Likewise, the analysis holds for fields of 
higher spin, where additional indices have been suppressed, as only the statistics of the 
field is relevant to this analysis. 

The operators obey the canonical commutation relations. 





cik) ^k' 


cr 


— 5k,k 




&k, bl> 


(J 


= l^k.k 




Sk, 6k/ 


(J 


= 0, 


where a = -|-1 for bosonic 


operat 


ors 


and a 



flkO-k/ 



(Tcij^/ak 



The free Hamiltonian is 
Ho 



k 



t^k 



(2.3) 
(2.4) 
(2.5) 



(2.6) 
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The iia,b are chemical potentials for particles and anti-particles respectively. If we have a 
U(l) symmetry, as in the case of a complex relativistic field, then = — A*6- For a real 
scalar field /Xa = /^fe = 0. The N^'^ are conserved number operators whose normalised 
eigenvectors are denoted by |nk) for particles, and |mk) for anti-particles. They are 
complete and satisfy 

iV^l^k) =nk|nk). (2.7) 

For systems described by Bose- Einstein statistics, nk = 0, 1, 2, . . ., whilst for Fermi-Dirac 
systems, rik = 0, 1. 

Thermal expectation values are denoted by angular brackets, and are defined by 
«...» = ^Trje-'^^o...} 

= ^ E ({nk},{mk}|e-^^°...|{nk},{mk}), (2.8) 

where ... means any operator, and P is the inverse temperature: P = 1/T {ks = 1). 
The trace has been taken over the number operator eigenstate basis [{jT-k}, {"^k}) = 
l^ki) ® ® ■ ■ ■ |?7T.ki) ® |?7?.k2) Zq\s the free partition function of the system: 

Z, = Tr{e-^^°} 

= (n^.")(n4)- (2-9) 

Z^^ is related to the Bose-Einstein and Fermi-Dirac distributions 

^a,b ^ j^g/3(c.k+Ma.6) _ (2.10) 

by 

Z^^^(l + N^^Y. (2.11) 

The following identities are useful for calculating thermal expectation values of sets 
of creation and annihilation operators 

00 

E nke-'^('"''+'^-.^)"'' = N^^ (1 + N^^) , (2.12) 

"k=0 

E {n],fe-^^^^+^''^'''>^ = N^"" (1 + K'') (1 + 2A^^'^) . (2.13) 

Thus for fields of either statistic, 

<aj,ap> = iV«4,p (2.14) 
<akaj,> = (1 + (jA^^) 5k,p (2.15) 

We now split the field operators arbitrarily into two parts 

V^(x) = '0+(x) + i)~{x). (2.16) 
5 



We call the parts "positive" and "negative" even though in general the split is not into 
positive and negative energy waves. In order to avoid possible confusion between f and 
+ superscripts, the hermitian conjugate of an operator is denoted hy ip = ip'^ . We do not 
assume that the operator and its hermitian conjugate are split in the same way. That is, 

we do not assume that ([V^]^) = This means care must be taken with the notation; 
we always have = [V'^]'*' rather than it being equal to [^/'^j^. 

To allow for products of more than one type of field operator, we now use the notation 



l/ji = 'IpiiXi). 



(2.17) 



Of course we will often be looking at products where some or all of the operators ipi{x) 
are equal or related by hermitian conjugation. 

Normal ordering is defined in terms of the arbitrary split ( p.l6|) , so generalising the 
traditional T = definition which is expressed in terms of annihilation and creation 
operators. We always define normal ordering to strictly mean that all + fields are placed 
to the right of — fields; otherwise the order of the fields is left unchanged. (So it differs 
from that used for fermions in [|10]-) For example, 

N [^i7/>2] := iVi,2 = ^ti^t + i'l^t + <^'4^2^t + A" ^2"- (2.18) 

Here ipi and ■?/'2 are any two fields with different space-time arguments (ti, xi) and (t2, ^2) 
respectively (which for convenience we usually suppress). The fields may be related, 
e.g. one often has ipi = ip{ti,:>ciy ,ip2 = ip(t2,^2)- For general splits A''i^2 may not be 
symmetric; Ni^2 7^ cr-^2,i- 

Time ordering is defined in the usual way so for the two-point case we have 



T := Ti,2 = ^1^2^(tl - t2) + (Tij2i^ie{t2 - tl). 



(2.19) 



It is straightforward to generalise this time ordering to path ordering as used in some 
finite temperature analysis. Then the theta functions of time in ( 2.19 ) become defined in 
terms of a path in complex time [Q]. As this does not effect the subsequent analysis we 
will not complicate our notation further. 

Wick's theorem is an operator identity. For two-point functions it relates Ti 2 to A^^i 2 
by defining the contraction Z}['0iV'2] = D12: 



Tl,2 



Nl.2 



o{t2 - ti) { ut, i^t] + k", V'2^1 + k", ^2~1 } • 



(2.20) 



(2.21) 



We will only consider splits, ( |2.16| ), for which D12 is a c-number - that is splits linear in 
the annihilation and creation operators. Note that the contraction, and hence the normal 
ordered product, is only symmetric, Di^2 = cr-D2,i, if 



^t,i^2 



+ 



i^i ,^2 



0. 



(2.22) 
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To higher order, Wick's theorem can be obtained from the generating functional 



T 



exp{-z / d^x ji{x)iJi{x)} 



1 



N 



exp{-z / d'^x ji{x)tpi{x)} 



X exp{-^ / d^x d% 3i{x)D[tlji{x)iljj{y)]jj{y)}, 



(2.23) 



where ji{x) are sources. Functional differentiation with respect to these sources gives 
Wick's theorem, which after using the symmetry properties of the time-ordered product 
is of the form 



2,...,2n 



E(-irE 



perm 



m=0 



a2,...,a2ri 



(2m) 



1 / " 1 

(n-m)\ [ . 2 -'^-^'-^^ 



(2.24) 



The sum takes Oi, . . . through all permutations of 1, ... , 2n, and p is the number of 
times one has to interchange pairs of fermions in moving from 1,. . . ,2n to ai, . . .a2n 
order. The sum over permutations can be simplified (even for non-symmetric normal 
ordered products and contractions) to give the usual form of Wick's theorem. For four 
bosonic fields, ( |2.24[ ) gives the well known result 



Ti 



2,3,4 



A^l,2,3,4 + /^l,2iV3,4 + /^l,3iV2,4 + /^l,4iV2,3 + /^2,3A^1,4 + 
^2,4A^1,3 + ^3,4A^1,2 + ^1,2^3,4 + ^1,3^2,4 + ^1,4^2,3- 



(2.25) 



To order 2n one has 



^1,2,. ..,2n — ^1,2,. ..,2n 

+ -Di,2A'"3^4,.,,,2n + all other terms with one contraction 

+ -Di,2-D3,4A''5 6,...,2n + all othcr terms with two contractions 



-|- D12D3 4 . . . D2n-i,2n + all othcr possible pairings. 



(2.26) 



Once again, we stress that (|2.23|) - (|2.26|) hold whatever the decomposition of the fields 
into positive and negative parts provided the contraction is a c-number. They hold 
whether or not the contraction is symmetric. 



3 Two-point functions 

In this section we prove that it is always possible to define a split -^(x) = il!~^{x) + ip'^i^x) 
such that -C A''i^2 ^= 0. This then ensures that Di^2 Ti^2 ^= o"-D2,i. 

A general linear decomposition of iplx) and iplx) into positive and negative parts is 

^+(x) = J2V^ ((1 - /k)ake-^^-^ + F^bie''-^) , (3.1) 

k 

^-(x) = ^ {hd^e-"'-^ + (1 - Fk)6i,e^'^-^) , (3.2) 

k 
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=Y.y^ (^7k4e*'" + (1 - G^)b^e-''^-^) , (3.3) 

k 

^-(x) = ^ ((1 - ^7k)4e''" + Gk6ke-^'=-") • (3.4) 

k 

Here = l/\/2uj^, and /k, (^k, and Fk are arbitrary complex numbers. We now 
impose the two independent conditions 

< N [tlJiiJ2] > = 0, (3.5) 

< N [V'2V'i] > = 0. (3.6) 



We find that (|3.5|) and (|3.6|) are satisfied if and only if 



/k^7k = -ctN^, 

(l-/k)(l-^?k) = (l + ^iV»), (3.7) 

FkGk = -crK, 

;i-Fk)(l-Gk) = (l + aiV^), (3.8) 



First consider ( |3.7] ) and ( |3.8| ) in the hmit T — 0. They are solved by /k = fi'k = Gk = 
Fk = which is just the traditional T = split as all the annihilation operators are in 

and the creation operators in ip~ . Taking the vacuum expectation value of ( p.20| ), 
the two-point version of Wick's theorem, we see that for T = 0, D12 = {0\T[ipiip2]\0) = 
Gi,2(T = 0), the zero temperature propagator. 

For T > 0, the equations for /k and (^k, -^k and respectively each have two solutions 

h = -aN^ + s''[K{K + a)]'/' , g^ = -aN^-s^[N^{N^ + a)Y/\ (3.9) 
F^ = -aN^ + s'[Nt{N^ + a)Y/' , = -aN^ - s''[Nt{N^ + (3.10) 

where s"''' = ±1 gives the two solutions in each case. Now, again, the thermal expecta- 
tion value of ( |2.2CI| ) gives D12 =<^ T['?/'i'?/'2] ^= ^1^2 (7"), where Gi^2iT) is the thermal 



propagator. 

For large k, we see that /k,5'k,-^k and Gk all tend to zero. Thus at high energies 
the split which gives zero normal ordered products at finite temperature is tending to 
be the same as the usual one used at zero temperatures. This is not suprising as the 
ultra-violet behaviour of finite temperature field theory is well known to be the same as 
zero temperature field theory. 

For bosons, ( pl9|) and ( p.9|) give real /k, g^, -Fk and G^, but ([V^]^) 7^ "ip^ unless T = 0. 
For fermions /k = and |/kp = Is'kP = ^k; ^^^1 similarly for Fk and G^. However, 
in this case, = ip~ for all T. 

For non-relativistic fields, we just set &k = 0. In the fermionic case, this recaptures 
the result of Thouless |TD|. 



For real relativistic scalar fields there is only one type of creation and annihilation 
operator. As Ok = &k and A^^ — ^k ~ -^k, a single Bose-Einstein distribution, one must 



have Fk = g^, Gk = /k in (|3.1| )- (|3^ ). This gives a general linear decomposition of a 
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real scalar field. Equations ( p.5| ) and ( p. 61 ) are now equivalent and their solutions are as 
above. 

Finally we consider the question of the symmetry of the normal ordered products for 
the case of two-point functions. On substituting (|3.1| ) to (|3.4| ) into ( |2.22|) one finds that 

Di,2 = ctD2,i ^ iVi,2 = ctN2,i ^f^ + g^ = 2f^g^ , F^ + G^ = 2F^G^. (3.11) 

This means that the split traditionally used in zero temperature work, where fk = gk = 
Fk = = 0, has a symmetric two-point normal ordered product. However (|3.9|) and 
( p.l(j| ) also satisfy ( ^.llj ). Thus the splits for which the thermal expectation value of 
two-point normal ordered products are zero, are guaranteed to have symmetric two-point 
normal ordered products. 



4 Four-point functions 

We now prove that the thermal expectation value of the four-point normal ordered prod- 
ucts A^^i, 2,3,4 vanish for the split of section 3. Hence we show that the four-point Green 
function decomposes into two-point Green functions for all V. For simplicity, we con- 
sider a real relativistic scalar field theory. Our results are generalised to other theories in 
section 5. 

Equation ( p.25|) can be used to relate iVi^2,3,4 to A^i.234; where in the primed normal 
ordered product the operators have been split according to the T = split {i.e. the one 
for which /k = fl'k = 0). The thermal expectation values of these two quantities satisfy 



< Nl,2,3A > = < Ar,,2,3,4 > - (^1,2 - Dl,2) (^3,4 " ^3,4^ 

Using the results of the previous section, it can be shown that 

D\ . - A,j = E(^k)^Ark \^e-'^-i^^-''') + e'*=-(^>-^^)] . (4.2) 

k 

The only terms which contribute to <^ 2 3 4 ^ must contain 2 positive and 2 negative 
parts. There are 6 such terms, and one of them is 

< ■^"(Xi)V^"(x2)V^+(x3)-^+(x4) > 

= i:k„k„k3,k4 V^,,V^,,V^,V^, (« ai^ai,&^,a^,. ») e^'=i-ie^^--^e-'=3-3e-.fc4 X4. (4.3) 

As in the introduction, we consider <C a|a2a3a4 3> where hi = St., and separate out 
explicitly the three contributions: 

<^ d|a2a3a4 ^ = 61^362^ -C a|a2aia2 ^ (1 — 61^2) 

+ 51,4^2,3 < 01020-20-1 > (1 - ^1,2) 

+ (5i,25i,3'53,4 < olaloioi > . (4.4) 
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Here 61^2 = <5ki,k2- From section 2, 



<C a[a2aia2 > 



A^iA^2 = <^ a[a2a2ai > 



(4.5) 



where iV, = A^ko and 



•C d|a|aidi ^ 



<^ d|aid|di ^ 
iVi (2iVi + 1) - iVi 

2(iVl)^ 



^ d|di > 



(4.6) 



Hence (^4.4|) reduces to 



(51,3^2,4 + 51,4^2,3) A^lA^2 

(5l,3^2,4 + 51,4(^2,3) < a{ai X 0,20,2 > 



(4.7) 



Given ( [4.7| ), ( |4.tj| ) and ( [4. 21 ) it is easy to evaluate the right hand side of ( [4.1D and find 
that 

< Ari,2,3,4 >= 0. (4.8) 

Thus taking the thermal expectation value of ( |2.25D , the four-point Green function 6*1,2,3,4 
is seen to decompose exactly into two-point Green functions Gi,2: 



Equation 



6*1,2,3,4 — G'i,2G'3,4 + Gi, 3(^2,4 + ^1,4(72, 3- (4-9) 

shows explicitly how contributions down by a factor of the volume, 
such as (|1.4| ), which is the last term of (|4.4| ), do in fact cancel. We just have to be careful 
to exclude the contribution of ( |1.4[ ) from ( |1.3| ) - hence the (1 — 5) terms in ( [4.4| ). It is clear 
from this derivation that there are no volume corrections to expectation value version of 
Wick's theorem 



5 n-point functions 



We now generalise ( ^.9| ) to n-point Green functions by using more powerful methods, 
namely the "trace theorem" and the path integral. Both of these methods may be applied 
to any type of field operator - bosonic or fermionic, relativistic or non-relativistic. The 
resulting generalised form of ( [4.9| ) then enables ( |4.8| ) to be generalised to all orders and to 
all types of fields. We also consider the symmetry properties of ra-point normal ordered 
products, and finally show how different splits may be related. 

5.1 The "Trace Theorem" 



This is described in detail in ||Il|, |12| and uses the cyclic property of the trace. 
Consider again 

>Vi,2,...,n =<^(a;i)?/J(a;2)...^(x„) >, (5.1) 
and write ^ 

^W = EE^kr/(k,x)d„(k). (5.2) 

k a=l 
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Here di(k) = and 62 (k) = aj^ (we still consider a real scalar relativistic field theory for 
notational simplicity). The thermal expectation value reduces to a trace over a operators, 
which we call S := Tr |e~'^^''da(ki)a5(k2) . . .at(k„)|. Since the commutation relations 
of the d's are known da(ki) can be commuted through d;,(k2), then ad^^,), and so on 
until it is finally commuted through d((k„). Then using the cyclic property of the trace, 
it can be taken to the left-hand side of the expression. Finally, after commuting through 
Q-0Ho Qj^g obtains 

l_e"»^-''i)5 = Tr{e-^^°d,(k3)...a<(k„)}[d,(ki),dfe(k2)] 

+ Tr {e-'3^«ab(k2) . . . at{K)] [«a(ki), de(k3)] 
+ ... 

+ Tr{e-^^«dfe(k2)...dt_i(k„„i)}[d„(ki),d4(k„)], (5.3) 



where Ai = —1, and A2 = 1. Now, one shows that 

5a,lh,2 < akiflka ^ +^a,2h,l < aj^i&ka > • (5.4) 



[da(ki),dfe(k2)] , , . -.t A ^;it 



Repeated use of ( [5.3| ) and ( p^ ) enable Wi,2,...,n to be expressed as a sum of products of 
the form Wi,2,...,n-2VVn-i,n- Iterating the argument, one obtains a generalisation of ( fl.9|) , 
that is 111 

(^^"(Xi, . . . , a;2„) = ^ f — j'-GaLaa • • • ■^G^a2n-l,a2n ) • (5-5) 
perm ^ ^ / 

The sum takes ai, . . . a2n through all permutations of 1, ... , 2n. Clearly this argument 
has no reference to volume. 



5.2 The Path Integral 



Equation (|5.5| ) may be obtained by using the path integral approach to thermal field 
theory. Let C be a directed path in the complex time plane whose end is —ij3 below the 
starting point, and let Tc denote path ordering of the operators with respect to C. The 
generating functional is then defined by 

Zo[j] :=Tr|e-^^°Tcexp||^f/t j rf^x j(t, x)0(t, x)|| , (5.6) 

where j(t,x) are classical sources coupled to the classical fields 0(t,x). This can be 
rewritten as 

Zo[j] = exp ^ dtdt' j d'lcd^x' Jit, x)Ac(a; - x')j{t', x')} . (5.7) 

Here Ac{x — x') = Gi^2 is the two-point thermal Green function or propagator - it is 
the Green function of the classical equation of motion for the fields 0(t, x) subject to the 
KMS boundary condition [^. As usual, the higher order Green functions are given by 



G"(xi,...,a;„) 



(5. 

j=0 



and so ( |5.5| ) is once again obtained. Again there are clearly no volume corrections of the 
type mentioned in the introduction. 
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5.3 Expectation values of normal ordered products 

We now consider thermal expectation values of normal ordered products of any type of 
field operators. Since both the path-integral and "trace theorem" can be generalised to 
all types of field operators B |TT|, ([5.5|) can be generalised to give 



G'-{x„ = E (-1)' (^^lGa,,a, • • • . (5.9) 

perm KU. Z Z J 



However, from ( |2.24| ) we know the relation between time-ordered products, normal or- 
dered products and contractions. Suppose we work with the splits satisfying (|3.71 ) and 
( p.8[ ) where the expectation value of the two-point normal ordered products is zero and 
the contraction is equal to the thermal propagator. Then (|5.9| ), with ( p.25|) for the four 
point case, imply that <^ A'^i, 2,3,4 ^= 0. This generalises the result proved directly for 
real relativistic scalar fields in section 4. 

Now we combine these results for two- and four-point functions with the six-point 
versions of (^]9|) and ( |2.24| ); they imply that -C -/Vi_2,3,4,5,6 3>= 0. Repeating, we see 
that the splits satisfying ( p. 71) and ( p.8|) ensure that the thermal expectation values of all 
normal ordered functions are zero, 

< iVi,2,...,2n >= 0. (5.10) 

All of these results are independent of volume. 

5.4 Symmetry of normal ordered products 

By definition the time-ordered function has the symmetry property Ti 2,...,n = (— 1)^7"^ 6 t 
where a,b, . . . ,t is any permutation of 1,2, ... ,n and p is the number of time pairs of 
fermionic fields are interchanged in moving from an order 1,2, ... ,n to an a,b, ... ,t order 
of fields. Using ( |3.11| ) it can be shown that 

(A,2 = aD2,i ^ N,,2 = ^iV2,i) iVi,2,...,„ = i-lYNa,t,...,t. (5.11) 
Thus for both the T = split and that of (|3.9|) and (|3.10D, the normal ordered products 



have the same symmetry as time ordered products. 



5.5 Relating two different splits 



Equation ( |2.23|) may be used to derive a relationship between different splits, i.e. between 



different types of normal ordering. Suppose N and A^' label two different kinds of normal 
ordering (for which D and D' are the corresponding c-number contractions). Then from 
( p.23| ) they are related by 



N' 



exp{-i / d^x ji{x)i)i{x)} 



N 



exp{-i / d'^x ji{x)i)i{x)} 



X exp{-i I d'xdS Ux) {D[^P,{x)^^{y)] - D'[^,(a;)^, (y)]) j,(i/)}. (5.12) 
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After functional differentiation tfiis yields 



perm \ )' 



perm 



n 1 



(5.13) 



The difference between normal ordering is therefore related to the difference between 
contractions. Now let D' and A^' refer to the usual T = split. Then D'12 = G'i,2(^ = 0) 
and -C A^' 0. Thus in calculating the thermal expectation value of a n-point time 
ordered function using ( p. 26 ), the terms containing <^ N' ^ are non- vanishing. (Recall 



that for this T = split normal ordered products are symmetric, and so the sum over 
permutations in ( |5.13|) simplifies greatly.) 

Now compare the T = split to the splits satisfying ( p.7| ) and ( |3.8| ). Let these 
splits have corresponding normal ordering and contraction labelled by A^ and D. Then 
•C A^ 3>= 0, and taking the thermal expectation value of ( |5.13| ) gives 



« Ar',,,...,,„ »= n(^-^')a.,-,.a., = E V^TT n(G(r)-G(r = Q)U_^,a2, 

perm ^ j=l perm ^ "" j=\ 

(5.14) 

Hence if one works with the usual T = split, the additional terms coming from ^ 
A^' are accounted for precisely by the thermal corrections G{T) — G{T = 0) to the 
zero temperature propagator. It may therefore make more sense to work with the splits 
defined by (^) and (pj). 

We can use the formula ( [5. 141 ) to provide a different view of several aspects of thermal 



field theory. For instance in real-time approaches there is a doubling of the degrees of 
freedom P, |, ||, ||, ^ 3- It is only the finite temperature corrections to the zero 
temperature propagators which mix the 'physical' and 'unphysical' degrees of freedom. 
Thus one could look on the doubling as being necessary only to take account of the normal 
ordered products, A^', where they are defined in the usual manner of zero temperature. 



6 Conclusion 

Contrary to claims in the literature, we have proved that for fields of arbitrary statistics 
there are no finite volume corrections to the usual relationship between thermal Green 
functions ( |5.9| ). Thus the operator approaches and path integral approaches to thermal 
field theory are indeed consistent in this respect. So in trying to understand the physics 
of the quark-gluon plasma formed in relativistic heavy-ion collisions, we have to look 
elsewhere to understand how finite volume affects the physics. 

In proving the above, we have defined and used a linear split of the field operators 
into positive and negative parts such that <^ A'i^2,...,2n ^= 0. This split enables usual 
T = field theory to be mimicked. By choosing to work with a split where such a result 
holds, we will simplify practical calculations. 

For instance when we are working with unphysical degrees of freedom at finite temper- 
ature, such as appear in gauge theories, one usually thermalises all the degrees of freedom. 
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However, it has been noted that such degrees of freedom need not be thermahsed . In 
the language of this paper, we would say that one may choose a different split (|2.16|) for 
each field. In particular, we may optimise the splits depending on how each field acts on 
the states in the ensemble. So with unphysical modes and an ensemble average taken over 
just physical states, the best split, normal ordering etc. for the unphysical modes is the 
zero temperature one, whereas we normal order the physical modes in terms of the usual 



thermal split of (pi]} and (|3.8|) . This idea can extended to out-of-equilibrium situations 
where different physical degrees of freedom may well be present in the ensemble in quite 
different amounts suggesting that we define normal ordering and contractions differently 
for each field so as to ensure that -C -/Vi,2,...,2n ^= always. 

We have refered only to the path-ordered approaches to thermal field theory, be they 
real- or imaginary-time formalisms. Another simple method of obtaining ( |5.9| ) is to use 
TFD 0, H, |, 1^, |[. There one must normal order with respect to the "thermal quasi- 
particle operators" C,,^,^^ and which annihilate the thermal vacuum, as this guarantees 
that all thermal expectation values of normal ordered products vanish, for example see 



19| . (For definition of the terminology and notation, see 0, [§.) Equation ( p.9|) then 



follows directly from ( |2.26| ) whereas in p] a rather lengthy derivation is given. This 
argument is also clearly independent of volume as it mimics the standard zero temperature 
field theory case. It is immediately obvious that these comments regarding TFD and 
normal ordered products apply to any situation where a Bogoliubov transformation is 
used. 
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